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Abstract 

We give a brief review of interacting NADE model in non-flat universe, we study the effect 
of spatial curvature 0,^, interaction coefficient a and the main parameter of NADE, n, On EoS 
parameter Wd and deceleration parameter q. We obtain a minimum value for n in both early and 
present time, in order to that our DE model crosses the phantom divide. Also in a closed universe, 
changing the sign of q is strongly dependent on a. It has been shown that the quantities Wd and 
q have a different treatment for various spatial curvature. At last, we calculate the statefinder 
diagnostic and w — w' analysis in non flat universe. In non flat universe, the statefinder trajectory 
is discriminated by both n and a. 
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I. INTRODUCTION 

The observational data of la supemova(SNIa) Jl|, large scale structure (LSS) |^ and 
cosmic microwave background (CMB) anisotropy [3] show that the universe is undergoing 
an accelerating expansion. It is believed that a dark energy with negative pressure can 
drive this accelerated expansion. The dark energy (DE) problem attracted a great deal of 
attention in the last decade. Up to now, many models such as ACDM , the models with 
a scalar field and modified gravity had been proposed l4|.Also in the last decade the other 

n n 

models base on quantum field theory such as holographic [5[ and agegraphic (ADE) [6| dark 
energy models are suggested. The latter is arisen from combining quantum mechanics with 
general relativity directly. It is worthwhile to mention that up to now, a completely successful 
quantum theory of gravity is not available. There are two main problems in dealing with 
ACDM model which are 'fine-tuning' and 'cosmic coincidence' problem j?]. The ADE model, 
which is proposed by Cai jo], is based on the line of quantum fluctuations of spacetime, the 
so-called Karolyhazy relation 6t = Xtp^^t^^^, and the energy-time Heisenberg uncertainty 
relation Est^ ~ t^^- These relations enable one to obtain an energy density of the metric 
quantum fluctuations of Minkowski spacetime as follows jst 

E5fi 1 ml 

Throughout this paper, we use the Plank unit [h = c = ks = 1) , where tp = Ip = l/rrip are 
Plank's time, length and mass, respectively. In ADE, this energy density is considered as 
density of dark energy component, p^, of spacetime. By considering a Friedmann-Robertson- 
Walker (FRW) universe, due to effect of curvature, one should introduce a numerical factor 
3^2 in (dJ |9|,Q. 

By making a model that considering DE, independent of the other mater fields, one 
can study the evolution of characteristics of dark energy of spacetime very well. The recent 
observational data from the Abell Cluster A586 supports the interaction between dark matter 
a„d dark energy Q. However, the strength of this „rteract,o,r ,s ,rot exactly .deirffied Q. 
Also, nowadays, many authors are interested to consider non-flat FRW universe [13|,lli. The 
tendency of a closed universe is shown in a suite CMB experiments [15i]. Besides of it, the 
measurements of the cubic correction to the luminosity-distance of supernova measurements 



reveal a closed universe 



16| . In accordance of all mentioned above, we prefer to consider a 
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model including dark matter and dark energy for a non-fiat FRW universe. 

The interacting new ADE (NADE) which is a new version of ADE with a conformal time 
instead of a cosmic time in FRW metric. The motivation of this new model is for that the 
original ADE model cannot explain the matter-dominated era jl?!] . Recently, the interacting 
NADE model has been investigated and some cosmological quantities such as deceleration 
parameter q, evolution behavior of fractional dark energy density fl'^, and equation of state 
(EOS) parameter wa = Pd/ Pd, are obtained isl. fiojl. 

The next step beyond Hubble parameter H = a/a and q = —d/aH"^, is to consider a 
new quantity contains a. A pair quantities which have been introduced by Sahni et al. and 
Alam et al. ||20i], are called statefinder pair {r,s}, as 

" 3(g-l/2)' ^ ^ 

The statefinder pair is a geometrical diagnostic tool which is constructed directly from a 
spacetime metric. The importance of such pair is to distinguish of the cosmological evolution 
behaviors of dark energy models with the same values of Hq and go at the present time. At 
future by combining the data of Supernova acceleration probe (SNAP) with statefinder 
diagnosis, we may choose the best model of dark energy. Up to now, many authors have 
investigated statefinder trajectories for standard ACDM model and quintessence 
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interacting quintessence models 
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221, 
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24| , chaplygin gas, the holographic dark energy models 



261], the holographic dark energy model in non-flat universe 27[, the phantom model 



28|, the tachyon 



the ADE model with and without interaction 



and the interacting 



NADE model in flat universe [19]. They had shown the stateflnder diagnosis is a useful 
tool for discrimination between various dark energy models. In addition to the stateflnder 
geometrical diagnostic, the another tool to distinguish between the different models of dark 
energy is w — w' analysis which is used extensively in the literature 
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30|. 



In this paper, in addition of cosmic behavior investigation, we study the stateflnder 
trajectories and w — w' analysis for interacting NADE model in a non-flat FRW universe. 
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II. COSMIC EVOLUTION IN NON-FLAT UNIVERSE WITH INTERACTING 
NADE 



As we mentioned in Sec. (jl]), the energy density can be defined in ADE model as 



Pd = (3) 
where the cosmic time T is defined as the age of the universe 

Introducing a conformal time rj which is defined as dt = adrj, the FRW universe is modified 

as 

ds'^ = dt^ — a^dx^ = a^^drf — dx^). (5) 

Therefore, by substituting the time scale r) in Eq.([3]), one can obtain the energy density of 
NADE model as 

_ _ f dt _ r da _ . _ 1 

^'^ 7f ' ^ J o, Jo Ha^' ^ a 



The corresponding fractional energy density is 

The Friedmann equation of a non-fiat FRW universe containing a new agegraphic dark 
energy and pressureless matter (baryons and dark matter) is 

where k = 1,0,-1 is curvature parameter corresponding to closed, fiat and open universe, 
respectively. Some recent observations reveal a closed universe with a present small fractional 
energy density Qko = ^/Hq ~ 0.02 [31]. Also we can write the Friedmann equation ([H]) in to 
another form with respect to fractional energy density Qi = Pi/pc, with pc = SmpH"^. Then 
we have 

nm + ^d = l + ^k- (9) 

The continuity equations including an interaction term between dark matter and dark energy 
become 

p'm + "iHpm = Q, (10) 
Pd + 3H{pa+Pd) = -Q, (11) 



where pd is dark energy pressure which is given by equation of state (EoS), Wd = Pd/pd- 

m, y, liol are 



Three forms of Q which have been extensively used in hteratures 



Q = ^aimlH^Vti- 



fid] « = 1 



i = 2 



m ) 

^d ~l~ ^mj ^ = 3 



(12) 



Differentiating Eq. ([7]) and using Eqs. ([8]) and ([9]), the derivative of Qd can be calculated 



as 



H 



H 



h - — - 

H'^ na 



na 2^ 2 2' 



(13) 
(14) 



where prime denotes the derivative with respect to In a and Qc = Q/Hpc = Saifli. The 
relations ( IT3ll and ( IT^ . also has been obtained in |l3| for third interaction form of Q. 
Using the relation ( !T4|) in ( !T3|) . we obtain a normal differential equation for Qd as 



f2 J — f^n 



/a 



;i_Q,)(3_2^) + nfc-Q, 
na 



where is given by 



a7(l - 



1 — 07 



Hear 7 is satisfied in the following equation 



— = a—- = 07. 



(15) 



(16) 



(17) 



From Eqs. (16]), ([T]) and ( ITTl) . the EoS parameter can be obtained as 

2v^ 



Wd 



-1 - Qcd + 



3na 



(18) 



where Qcd = Qc/(3fid). The evolution behavior of Wd (fTSj) by using (fTSl) . is given by 



+ 



3r;,a 



1 + fifc - 3(]d - 2^-^(1 - Qd) - Qc 
na 



-Q'c 



cd' 



(19) 



Also the total EoS parameter is obtained as 

_ ptot _ ^dWd 

Wtot — — , , 

Ptot -L + "fc 



(20) 
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To achieve an accelerated expansion, it is required that Wtot < —1/3. Therefore at present 
time in a closed universe (a = l,Qko = 0.02, Qdo = 0.72), the minimum value of n with 
Q!j = 0.1 for various forms of Q, can be obtained as 



n> < 



0.901; i = 1 
0.993; i = 2 
0.845; i = 3 



(21) 



In interacting NADE, based on Eq. fllSp , we see that at present time, wa can cross the 
phantom divide {w < —1) provided 



na > < 



0.566 
1.358 
0.399 



1 
2 

3 



(22) 



The third value of fl221) . has been also expressed in 13|. The relation fl2^ show that for the 
same value of a, Ui is increasing by moving from i = 3 in to z = 1 and then i = 2. 



The deceleration parameter q can be calculated from Eqs. (fT3|) and (fT5|) as 



na 



2 2 



Qc 

2(1 - 07) 2 



(23) 
(24) 



III. STATEFINDER DIAGNOSTIC OF NADE IN INTERACTING NON-FLAT 
UNIVERSE 

Now we find the statefinder pair {r,s}, which was expressed in SeclH From the definition 
of q and H, the parameter r ([2]) can be written as 

r = ^,-3q-2. (25) 

Using Eqs. ([HD, ([18]) and ([23]) we have 

^ = 1 + ^QdWdiwd + Qcd + 2) - In^w'^ + ^Qk- (26) 
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Hence, the Eq. (125|) can be obtained as 



9 3 
r = 1 + fifc + -^^dWdiwd + Qcd + 1) - -^^dw'^. 



(27) 



In a non flat universe, Evans et al. 



32| generalize the definition of parameter s ([2]) as 
_ r - Qtot 



(28) 



where the total fractional energy density is VLtot = + = 1 + ^k- Therefore from this 
new definition we have 



S = I + Wd + Q 



cd 



(29) 



The relations (I26D,(I27D and ([29]) reduced to (35), (37) and (38) of Ref. |l9| in the limiting 
case of flat universe. By omitting w'^ between ( 127|) and (|29|) . we find r in terms of s as follows 
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r = 1 + fifc + -sildWd- 



(30) 



IV. NUMERICAL RESULTS 



In this section, first we give the complete numerical description of the NADE model and 
then examine the NADE model with statefinder diagnostic tool and w — w' analysis. Here 
we consider the first case of interaction form, Q = Qi = dai'm'^H^fld, in ffT^ with a = ai. 
In this case, the differential equation for Qd (II5D can be reduced as 



n', = Qdlil - ^d){S - ^^) - SQda + fifcj 



na 



(31) 



Substituting Qi in Eqs. (fT8|) . ( [T9|) . (!27|) and (129|) . yields the following equations: 



, , 2y^d 

Wd = -1 + — a 

6na 



[(1 + (^,.)- 31^.(1 + a) -^fl 



"ina 



na 



9 3 
r = l + VLk + -^dWd{wd + a + 1) - -^dw'^- 



I + Wd + a - 



Wd 

3wj 



(32) 

(33) 
(34) 
(35) 



Because of vanishing fi^ at a — )■ 0, the second term of Eq. (l32|) can be larger or smaller than 
unity at this time. This term is examined by numerical analysis and we see that it is much 
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lower than one for higher values of model parameter n. Hence, NADE model can cross 
the phantom divide {wd < — 1) at a — t- 0, independent of the contribution of the spatial 
curvature of the universe. 

On the other hand, at the late time (e.g., a — t- oo) when fi^^ — )■ 1, tends to —1 and 
NADE model mimics the cosmological constant. Here we focus on Eg. fl32p in more details. 
Wd crosses the phantom divide {wd < — 1) when 2y/rLd/'ina < a, otherwise Wd can not cross 
the phantom divide. In SectUll we obtained a condition for na > 0.566 at present time for 
crossing the phantom divide. For the best value of n = 2.7 which has been obtained from 
astronomical data for NADE model, we should set a > 0.2 to have < — 1 at present time. 
However, the situation is changed at the early time, since both a and Qd tend to zero at that 
time. For example the values n = 4 and a = 0.1 satisfies the relation 2^/^ld/3na < a at the 
early times and Wd crosses the phantom divide, while the condition 2^/^ld/3na < a can not 
be satisfied for the values n = 3,a = 0.1 and the phantom divide can not be achieved at the 
early times in this case. 

Fig.(l) shows the evolution of Wd of interacting NADE model in terms of scale factor. In 
Fig.(l-a) we illustrate the evolution of Wd for fixed model parameters n = 4 and a = 0.1, 
in open, flat and closed universe. All three cases give the phantom divide at the early time 
and mimic the cosmological constant at the late time. In Figs.(l-b,l-c), the dependence 
of the evolution of Wd on the model parameters n and a are investigated. Here we choose 
the closed universe with the present spatial curvature Qko = 0.02. In Fig.(l-b), by fixing 
a = 0.1, we vary the parameter n as 3, 4 and 5. In the case of n = 3 the interacting NADE 
model can not cross the phantom divide, while for n = 4 and n = 5 the phantom divide 
can be achieved. In Fig.(l-c) we fix n = 4 and vary a as a = 0.0 , a = 0.1 and eventually 
a = 0.15. It can be seen that the phantom divide is archived for a = 0.1 and a = 0.15 and 
it can not be access for a = 0.0 (the NADE model without interaction). 

The other cosmological parameter which we demonstrate, is the deceleration parameter q. 
The parameter q in NADE model for non-flat universe is given by Eq.( |23ll . In the early time, 
where fi^ — )• and Qk — ^ 0, the parameter q converges to 1/2, whereas the universe has been 
dominated by dark matter. In Fig(2), we show the evolution of g as a function of cosmic 
scale factor for different model parameters of NADE model and also for various contribution 
of spatial curvature of the universe. In Fig.(2-a), the dependence of the evolution of q on 
the spatial curvature of the universe {k = 0, 1, —1) is sketched for n = 4 and a = 0.1. 
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In this model, the deceleration parameter q crosses the boundary g = from g > to 
g < 0. This implies that the universe undergoes decelerated expansion at the early time 
and later starts accelerated expansion. The transition from decelerated expansion to the 
accelerated expansion occurs gradually from closed, fiat and open universe. However, the 
difference between them is very little, but we can interpret that the transition occurs earlier 
in closed universe. In Figs.(2-b,2-c) the evolution of q in terms of the scale factor is plotted 
for different values of n and a in the case of closed universe with Vlko = 0.02. In Fig.(2-b), 
we set a = 0.1 and change n as 3, 4 and 5. Here the change on the sign of q is taken place 
at similar a for all values of n. It should be noted that in the accelerated universe (g < 0), 
the parameter g is smaller for higher values of n while in the decelerated universe (g > 0), 
g is larger for higher values of n. In Fig.(2-c), by fixing n = 4, and changing a as 0.0, 0.1 
and 0.15 the behavior of g is studied. Here the universe starts accelerated expansion earlier 
when a is more. 

At following, we calculate the evolution trajectories in the statefinder planes and analyze 
the interacting NADE model in non-fiat universe with statefinder point of view. The stan- 
dard ACDM model in non fiat universe corresponds to the fixed point (s = 0,r = 1 + Qk) 
in the s — r plane [32]. One way to test the ability of a given dark; energy model is the 
deviation value of the model from the fixed point (s = 0,r = 1 + Qk) in s — r diagram. Let 
us start with Eqs.f l33|l35] and which describe the evolution of statefinder parameters r, 
s and also w'^. It is easy to see that in the early time, w'^ — )■ 0, and from ( [341) and (!35ll . 
s — )■ and r — )■ 1. Also it is worth to estimate the values of u'^, s and r at late time when 
a — 7- oo. From Eqs.([3lj) and(l35l). we obtain s ~ a and r ~ 1 at late time. So the statefinder 
parameters (s,r) reach to the (a,l) at the late time. 

In Fig. (3) the evolution trajectories of statefinder of interacting NADE model with inter- 
action form of Qi is plotted. In Fig.(3-a) , the evolution trajectories is plotted for different 
closed, fiat and open universe for fixed parameters n = 4 and a = 0.1. The present values 
of the statefinder parameters sq and tq is denoted by circle {sq = 0.154,ro = 0.535), star 
(so = 0.154,ro = 0.515) and square (sq = 0.153,ro = 0.497) symbols for closed, fiat and open 
universe, respectively. It should be noted that the evolution trajectories start form fixed 
point (s = 0,r = 1) at the early time, as mentioned above. It can be seen that the different 
curvatures will lead to different evolutionary behavior in the statefinder plane, starting from 
the same fixed point {s = 0,r = 1) at the early time. The curvature will affect the today's 
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value of statefinder parameter. Fig.(3-a) shows that the distance to ACDM fixed point in 
closed universe is shorter than of obtained in flat universe and both of them is shorter than 
that distance in open universe. Also in closed universe the value of r is the largest, while s 
is the smallest at present time. 

Figs.(3-b & 3-c) indicate the dependence of the evolution trajectories of stateflnder di- 
agnostic on n and a in closed universe. Fig.(3-b) shows the influence of the variation of 
n on the evolution trajectories of stateflnder for flxed parameter a = 0.1. Here the evo- 
lution trajectories is calculated for n = 3,4 and 5. Symbols on the curves represent the 
present value of stateflnder. The circle indicates (sq = 0.205,ro = 0.406) for n = 3, the 
star symbol indicates (sq = 0.154,ro = 0.535) in the case n = 4 and the square denotes 
{sq = 0.123,ro = 0.620) for n = 5. Increasing the parameter n will lead to shorter distance 
between present values (scTq) and flxed ACDM in this diagram. Also we can see that the 
higher value of n makes the larger value of r and smaller value of s. 

In Fig. (3-c), we redo the previous calculation in Fig.(3-b) for flxed parameter n = 4 
and varying parameter a as 0.0, 0.1, 0.15. The variation of a also change the evolu- 
tion trajectories of stateflnder. The present values of stateflnder for different values of 
a are: [sq = 0.162,ro = 0.564) for a = 0.0, [sq = 0.154 ,ro = 0.535) for a = 0.1 and 
{sq = 0.151,ro = 0.521) for a = 0.15. We can see that the higher value of a makes the 
smaller value of r and also the smaller value of s. It is worth noting that in the case of flat 
universe, only the parameter a can discriminate the evolution trajectory in s — r plane and 
the parameter n can only sperate toady's value of s and r(see Figs. (3 & 4) of ISj). Here in 
non flat case, both parameters n and a can discriminate the evolution trajectories in s — r 
plane (see Figs.3-b & 3-c). 

At last, we study the interacting NADE model in non flat universe using the w — w' analysis. 
In this analysis, the standard ACDM model corresponds to the flxed point {wd = — l,w^ = 0) 
in the w — w' plane. The evolution of Wd and w'^ is given by Eqs.(!32l [33l) . In Fig. (4) the 
evolution trajectories of in w — w' plane is shown for different parameters and various 
curvatures. In Fig.(4-a), we show the evolution trajectories of w'^ , by flxing the parameters 
n = 4 and a = 0.1 for different spatial curvatures. Here we see that the various spatial 
curvatures gives the different evolutionary behavior in the w — w' plane. The evolution 
trajectory of different curvatures converges to the flxed point {wd = —l.l,w'^ = 0.0) at 
late time. The present value of Wd and w'^ are: (w° = — 0.957, w'^ = —0.104) in closed. 
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(ujO = -0.957,w'^ = -0.106) in flat and = -0.957,^2' = -0.107) in open universe. 
In Figs.(4-b, 4-c), the evolution trajectories inw — w' plane is discussed for the case of closed 
universe. In Fig.(4-b), by fixing a = 0.1, we vary the parameter n as 3, 4 and 5. The today's 
value of w^, w'^ is denoted by symbols on the lines. The circle symbol shows the present 
value = -0.910, w^f = -0.142) in the case of n = 3, star (w^ = -0.957, w^J* = -0.104) 
for n = 4 and square {w^ = — 0.986, = — 0.082) for n = 5. Increasing the parameter 
n will lead to closing the present values of w^, w'^ tend to the fixed point (—1, 0) in w — w' 
plane. The higher value of n obtains the smaller values of Wd and the bigger value of w'. 
In Fig. (4-c), we fix = 4 and vary a as 0.0,0.1 and 0.15. The present values of 
w^,w'^ are shown by symbols on the lines. The circle symbol indicates the present value 
[w^^ = -0.857, w'^ = 0.051) in the case of a = 0.0, the star shows (u;^ = -0.959, w^^^ = 0.035) 
for a = 0.1 and the square represents {w'^ = —1.01,w'^ = 0.027) for a = 0.15. Increasing 
the interaction parameter a would lead to decreasing the present value {w^, w'^) in w — w' 
plane. 

V. CONCLUSION 

In this work, the interacting NADE model in non-fiat universe has been given. We studied 
the effect of spatial curvature flk, interaction coefficient a and the main parameter of NADE, 
n. On EoS parameter Wd and deceleration parameter q. We showed that in the early and 
present time for ai = 0.1, the phantom divide is not available for n = 3 and it is achieved for 
n > A. By increasing n and a in a closed universe, the trend of Wd decreases. We obtained 
a minimum value for n in both early and present time, in order to that the NADE model 
crosses the phantom divide. It was shown that the treatment of both parameter Wd and 
q are dependent on the type of spatial curvature. At last, we investigated the interacting 
NADE model in a non-fiat universe by means of statefinder diagnostic and w — w' analysis 
viewpoints. Here we showed that the spatial curvature can affect the evolution trajectories 
in (s,r) and {wd,w'^) planes. Also the trajectories in these planes can be affected by the 
model parameters of interacting NADE, n and a. In non fiat universe, the statefinder 
trajectory is discriminated by both n and a. It should be noted that in the case of flat 
universe, only the parameter a can discriminate the evolution trajectory in s — r plane and 
n can only discriminate toady's value of s and r(see Figs. (3 & 4) of \)^)- While in non flat 
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universe, both parameteres n and a in addition to discrimination of {s,r} at present time, 
can discriminate the evolution trajectory in s — r plane (see Figs.3-b & 3-c). It is worthwhile 
to mention that all computations is reduced to previous work in the limiting case of flat 
universe llQll. 




-0.4 
-0.5 
-0.6 

-0.7 
-0.8 
-0.9 
-1 
-1.1 



Fig.(l-b) 



Fig.(l-c) 





£1°=0.02 

k 


1 
1 
1 


a=0.1 


\ 

■ \ 








FIG. 1: The evolution of EoS parameter, Wd, versus of a for different model parameters n, a and 
different curvatures Vt^- Fig.(l-a): Dependence of Wd on , for n = 4 and a = 0.1. The dashed, 
solid and dotted-dashed lines represent the closed, flat and open universe, respectively. Fig.(l-b): 
the evolution of Wd for n = 3 (dashed line), n = A (solid line) and n = 5 (dotted-dashed line), 
by hxing a = 0.1, in a closed universe . Fig.(l-c): the evolution of Wd for a = 0.0 (dashed line), 
a = 0.1 (solid line) and a = 0.15 (dotted-dashed line) , by fixing n = 4, in a closed universe. 
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Fig.(2-a) 

0.6 1 




FIG. 2: The evolution of q versus of a for different parameters n and a and various contribution 
of spatial curvatures ^k- Fig.(2-a): Dependence of q on Q,k for values n = 4 and a = 0.1. The 
dashed, solid and dotted-dashed lines represent the open, flat and closed universe, respectively. 
Fig.(2-b): the evolution of q versus a, by fixing a = 0.1, for n = 3 (dotted-dashed line), n = 4 
(solid line) and n = 5 (dashed line) in closed universe. Fig.(2-c): the evolution of q versus a, by 
fixing n = 4, for a = 0.0 (dotted-dashed line), a = 0.1 (solid line) and a = 0.15 (dashed line) in 
closed universe. 
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FIG. 3: Fig.(3-a): Evolution trajectories of the statefinder in the r — s plane for open (dotted- 
dashed line), fiat (solid line) and closed universe (dashed line). Square, circle and star symbols on 
the curves are the today's values of the statefinder parameters {sq^tq). Here the model parameters 
are n = 4 and a = 0.1. 

Fig.(3-b): Evolution trajectories of the statefinder in the r — s plane for different values of n as 
3 (dottcd-dashed line), 4 (solid line) and 5 (dashed line) , by fixing a = 0.1 in closed universe. 
Fig.(3-c): Evolution trajectories of the statefinder in the r — s plane for different values of a as 0.0 
(dotted-dashed line) ,0.1 (solid line) and 0.15 (dashed line) , by fixing n = 4 in closed universe. 
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Fig.(4-a) 
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FIG. 4: Fig.(4-a): the evolution trajectories in Wd-,w'^ plane, by fixing n = 4 and a = 0.1, for 
closed (dashed line), flat (solid line) and open universe (dotted-dashed line). The present value of 
indicated by symbols on the curves. Fig.(4-b): Evolution trajectories in Wd,w'^ plane for 
different values of n as 3 (dotted-dashed line), 4 (solid line) and 5 (dashed line) , by fixing a = 0.1 
in closed universe. Fig.(4-c): Evolution trajectories in Wd,w'^ plane for different values of a as 0.0 
(dotted-dashed line), 0.1 (solid line) and 0.15 (dashed line) , by fixing n = 4 in closed universe. 
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